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We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as 
calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) 
type charge partitioning scheme accounts for intra-molecular contributions and overcomes the sys¬ 
tematic underestimation of the non-adiabatic losses obtained within the prevalent independent atom 
approximation. The quantitative agreement obtained with theoretical and experimental benchmark 
data suggests the LDFA-AIM as an efficient and reliable approach to account for electronic dissipa¬ 
tion in ab initio molecular dynamics simulations of surface chemical reactions. 
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A central challenge in energy and catalysis applica¬ 
tions is to transfer energy specifically into those degrees 
of freedom that actually drive a desired surface chemi¬ 
cal reaction - and to keep this energy in these degrees 
of freedom for a sufficiently long time. In this transfer 
of energy, losses due to electronic non-adiabaticity can 
be an important dissipation channel mm- In aiming to 
assess this channel for systems of technological interest, 
predictive-quality calculations would be a valuable addi¬ 
tion to experimental endeavors. Especially for chemical 
reactions at frequently employed metal substrates, how¬ 
ever, a corresponding methodology has not yet been well 
established. 

To date, most accurate solutions of the full nuclear- 
electron wave function are restricted to systems of the 
complexity level of gas-phase Eb" 1 " [3j. In the limit of 
weak non-adiabaticity as pertinent to electron-hole (eh) 
pair excitations during adsorbate dynamics on metal sur¬ 
faces, less rigorous approaches rely on mixed quantum- 
classical dynamics. The imposed computational burden 
nevertheless still restricts their practical use to simple 
metals and sub-picosecond time scales II 0, to sym¬ 
metric adsorbate trajectories 0 0, or to only qualita¬ 
tive accounts of the metal electronic structure 0 0- 
Presently, it is thus only the concept of electronic friction 
[HM3] and its incorporation into classical molecular dy¬ 
namics (MD) simulations on the Born-Oppenheimer po¬ 
tential energy surface (PES) Vpes [TH - ET] that promises 
predictive-quality and material-specific trajectory calcu¬ 
lations over an extended period of time. 

Particularly the local-density friction approximation 
(LDFA) [16) 22j and for molecular adsorbates an ad¬ 
ditional independent atom approximation (IAA) [[161 - 


m provide a further decrease in computational cost. 
This has allowed for first accounts of electronic non- 
adiabaticity in large-scale MD simulations based on a 
first-principles and high-dimensional description of the 
underlying PES - either interpolated [16H19I or even 
evaluated on-the-fly within ab initio MD simulations 
Ennui. However, due to the drastic simplifications in¬ 
troduced with the IAA, the validity of the LDFA for¬ 
malism for molecular adsorbates per se has been con¬ 
troversially discussed [16] [23] [24]. By construction the 
IAA does not resolve the electronic structure of the in¬ 
teracting molecule-surface system and in particular the 
location of the molecular frontier orbitals in the surface 
band structure. It can thus, for instance, not reproduce 
the enhancement of friction coefficients close to the tran¬ 
sition state of a molecular dissociation event on metal 
surfaces [25] . On the other hand, the necessity of an 
accurate description of such regions of enhanced friction 
for the overall non-adiabatic energy dissipation has been 
questioned, as the typically low velocities in these regions 
effectively suppress the contribution of the friction term 
within the dynamics 01- 

Despite the success in recent applications, it thus re¬ 
mains elusive to which extent the limitations of the preva¬ 
lent IAA carry over to actual observables. In this situ¬ 
ation, the vibrational lifetimes of high-frequency adsor¬ 
bate modes can provide a sensitive measure, as they are 
largely governed by energy dissipation in the electronic 
non-adiabatic channel E5H2TL Accurate experimental 
reference data then allows for a substantiated assessment 
of the quality of the non-adiabatic description. In this 
study we perform such an assessment, primarily focus¬ 
ing on the internal stretch mode of two systems which 
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have been studied most extensively and conclusively by 
experiments: CO adsorbed on Cu(100) and Pt(lll) [26] . 
Despite the largely different surface frontier orbital lo¬ 
cations and concomitant hybridizations at the transition 
and noble metal surface, we find the LDFA-IAA to al¬ 
ready exhibit a good qualitative performance with re¬ 
spect to experimental [23, J25] and theoretical m no 
benchmark data. Rather than an explicit account of 
the surface band structure, our analysis suggests missing 
intra-molecular contributions as reason for the remaining 
differences. Approximately incorporating such contribu¬ 
tions through a numerically efficient atoms-in-molecules 
charge partitioning indeed yields consistent lifetimes for 
a range of diatomic adsorbate systems. 

In friction theory, all non-adiabatic effects due to the 
excitation of eh pairs in the metal substrate are con¬ 
densed into a single velocity-dependent dissipative force 
that augments the classical equations of motion 


rrii 


d 2 R la 
d t 2 


dVpES 

dR ia 


N 3 

(i) 


Here, small latin and greek subscripts denote atoms and 
Cartesian degrees of freedom, respectively, and N is the 
total number of atoms of mass rrq and position R; in 
the system. The fluctuating white noise force JF, q (T) 
becomes negligible at very low temperatures and vanishes 
exactly at OK. Every element of the friction tensor rj is, 
in principle, a function of all nuclear coordinates. Within 
the spirit of weak coupling, the focus is usually on the 
diagonal contributions describing the electronic friction 
felt by each atom muz]. 

These atomic friction coefficients r]i a jp = 7]i a Sij6 a p can 
e.g. be calculated within the quasi-static regime building 
on time-dependent density-functional theory (DFT) as 
suggested by Persson and Hellsing mm- While insight¬ 
ful and generally in good agreement with experiments, 
the accurate numerical evaluation of this approach is 
challenging in practice and has hitherto been limited 
to low-dimensional potentials describing the adsorbate- 
metal interaction [Him sues EH. This shifts interest 
to more effective schemes and there in particular to the 
local-density friction approximation. The LDFA intro¬ 
duces isotropic scalar atomic friction coefficients 7? FDFA , 
such that ?7 ^ fa = r]f DFA SijSai 3 - These coefficients can 
be calculated very efficiently from the scattering proper¬ 
ties of an atomic impurity embedded in a free electron 
gas (FEG) [lOUniES]. In order to relate this model to 
the actual motion of adsorbates, the embedding density 
of the FEG p em b is then chosen as the electron density 
value of the clean metal surface p sur f at the position of 
the adsorbate atoms. In its application to molecular ad¬ 
sorbates this implies an independent atom approxima¬ 
tion which has been regarded and discussed as being 
inherently included in the LDFA formalism in this con¬ 
text [M119EI]- As a result of the IAA, the employed 


atomic friction coefficients are insensitive to the molec¬ 
ular nature of the adsorbate, yet can be evaluated very 
efficiently. The only input variable is the clean surface 
electronic density, which, assuming a frozen surface, has 
to be calculated only once in advance. 

We obtain the two-dimensional PES VpEs(dcOi Zoom) 
of an adsorbed CO molecule as a function of its bond 
length dco and center-of-mass (COM) height Zqom 
above the frozen surface through DFT calculations. Each 
PES is supported by 442 data points that are calcu¬ 
lated with the CASTEP plane-wave pseudopotential code 
[SB] and subsequently interpolated using bivariate cu¬ 
bic splines. Electronic exchange and correlation (xc) is 
treated on the GGA level in terms of the PBE func¬ 
tional m- The metal surfaces are modeled by five 
layer slabs with a separating 20 A vacuum distance. We 
consider top-site adsorbed CO molecules on one side of 
the slab within a c(2 x 2) and (v/3 x v / 3)R30° surface 
unit-cell on Cu(100) and Pt(lll), respectively. In both 
cases, the molecular axis is perpendicular to the sur¬ 
face, with the C atom coordinated to the metal atom. 
At the employed computational settings (600 eV cut-off 
energy, ultrasoft pseudopotentials [35], (10 x 10 x 1) 
and (11 x 11 x 1) Monkhorst-Pack k-point grids [39] for 
Cu(100) and Pt(lll), respectively) the PES data points 
are converged to < 5meV. Our investigation is not af¬ 
fected by the well-known CO adsorption puzzle and the 
concomitant wrong absolute depth of the adsorption well 
[40] . This is confirmed by essentially identical lifetimes 
we obtain when using PESs generated with a van der 
Waals-xc-functional [2] that leads to a stabilization of 
the top site [42] , 

Vibrational lifetimes r are extracted from classical 
MD simulations on the interpolated PESs by numerically 
solving Eq. Q. Within the LDFA the atomic friction co¬ 
efficients r]f“ A are calculated from the scattering phase 
shifts of the Kohn-Sham orbitals at the Fermi momen¬ 
tum Sf = Si(k f) for an atomic impurity embedded in a 
FEG of density p e mb,i uni uses] 

A 00 

^ LDFA (Pemb,i) = 7Tp ^ mh -‘ J2(l+ l)sin 2 [<5f - c5f +1 ] , (2) 
F 1=0 


where p* AA j = /9 sur f(R;) within the IAA as described be¬ 
fore. Assuming a constant energy dissipation rate and 
thus an exponential decay of the vibrational energy E v ib, 
the lifetime r can be extracted from the simulations by 
a logarithmic fit of E v ib versus time t. To initialize our 
simulations, we assign the adsorbate stretch-mode a pro¬ 
jected kinetic energy of Hu, where u is the normal mode 
frequency. Higher initial kinetic energies up to 5 Hu result 
in minute lifetime changes of less than 0.1 ps. 

Figure [l] shows the vibrational lifetimes that result 
from our simulations. We compare them to experimental 
values obtained from pump-probe spectroscopy [2E1 EHj 
and theoretical values as published by Forsblom and 
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FIG. 1. (Color online) Vibrational lifetimes for CO on (a) 
Cu(100) and (b) Pt(lll). Values as obtained within the 
independent atom approximation (IAA) and the atoms-in- 
molecules (AIM) approach from Eq. are contrasted to 
corresponding predicted lifetimes published by Forsblom and 
Persson (FP) [3Q: and Krishna and Tully (KT) [31]. For com¬ 
parison, experimental values as obtained from pump-probe 
spectroscopy by Morin et al. for CO on Cu(100) [29] and 
Beckerle et al. for CO on Pt(lll) [28] are shown as a dotted 
line and a blue stripe further indicating the reported experi¬ 
mental uncertainty. a CO on Cu(lll) 


Persson (FP) [3D], and Krishna and Tully (KT) 31]. 
The latter two are both based on the orbital-dependent 
Persson/Hellsing expression mentioned above [IT] [32], a ^" 
beit obtained from different derivations and relying on 
slightly different numerical treatment. For both systems 
the LDFA-IAA lifetimes agree fairly well with the theo¬ 
retical values from FP and KT, as well as with experi¬ 
ment. With all numbers falling within one order of mag¬ 
nitude, the current data thus does not further support 
the harsh criticism the LDFA-IAA was faced with be¬ 
fore [23]. Generally, they are instead consistent with the 
good LDFA-IAA performance reported in earlier studies 
on non-adiabatic energy losses of various ions scattered 
off metal surfaces [33] and on the vibrational damping 
of atoms on metal surfaces mmm- Conspicuously, 
however, in these studies on adsorbate atoms the agree¬ 
ment was even more quantitative and lacked the system¬ 
atic underestimation (overestimation) of LDFA-IAA non- 
adiabatic energy losses (lifetimes) apparent in Fig. [TJ 

Rather than from a generally insufficient account of the 
electronic structure of the interacting adsorbate-surface 
system, this suggests that LDFA-IAA deficiencies arise 
particularly in the treatment of molecular adsorbates as 
isolated adatoms. In this respect - within the underly¬ 
ing atomic embedding model - systematic shortcomings 
lie in the complete neglect of both adsorbate-substrate 
as well as intra-molecular contributions to p em b- Consid¬ 
ering only the density of the clean surface, p em b is sys¬ 
tematically underestimated; consistent with the under¬ 
estimated energy losses observed in Fig. [T] Even more, 
by construction a thus defined p em b is also blind to dy¬ 
namical changes of intra-molecular bond distances and 



FIG. 2. (Color online) (a) Change of the C-O bond length 
Adco on Cu(100) over the first few periods of the stretch 
mode, (b) Embedding density p em b,c of the C atom repre¬ 
sented by the Wigner-Seitz radius r s = (3/(47rp e mb,c)) 1 ^ 3 in 
atomic units (a.u.). (c) Corresponding carbon friction coeffi¬ 
cient. Solid red lines refer to the IAA, dotted orange lines to 
an AIM embedding density according to Eq. (j3|. 

characters. Figure [2] illustrates this for the vibrational 
motion of CO at Cu(100). The LDFA-IAA friction co¬ 
efficient of the carbon atom decreases with smaller C-0 
bond lengths dco and thus larger distances Zq of the 
C atom from the metal surface. This is due to the fact 
that the IAA only accounts for changes of p sur f along the 
vibrational coordinate. The LDFA-IAA thus effectively 
decomposes the molecular vibration into two indepen¬ 
dent adatom vibrations. Yet, even without any intra¬ 
molecular bonding there will be an increasing overlap of 
the O and C atomic densities with decreasing dco- This 
intra-molecular contribution to p em b at the position of 
every constituent atom is thus completely missed in the 
LDFA-IAA approximation. 

A straightforward way to account for such contribu¬ 
tions is to perform a charge decomposition through a 
projection scheme like Hirshfeld’s analysis [35] ■ For any 
given adsorbate configuration this provides the projected 
density of any adsorbate atom pf, and corresponding 
sharing function w] 1 at any position R [46] . Within an 
atoms-in-molecules (AIM) picture we can then define as 
embedding density for atom i at position R ; 

P^Sb,i = PSCF (R-i) - pf (R-i) 

= [l - wf (Rj)] Pscf(Ri) ■ (3) 

We thus consider as embedding density the full self- 
consistently calculated density Pscf of the entire inter¬ 
acting adsorbate-surface system just without the contri- 
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bution associated with atom i. This naturally contains 
density contributions from all other atoms in the system 
(both substrate and adsorbate) and consequently car¬ 
ries an implicit dependence Figure [2Tb) 

contrasts the correspondingly obtained to 

during the CO vibrational motion. Intriguingly, the two 
quantities differ not only quantitatively, but even exhibit 
a reversed phase behavior. The account of the density 
contribution of the O atom in p A ™c thus outweighs the 
influence of the clean surface density seen by the IAA: At 
the largest vertical heights of the carbon atom Zq , where 
PembC * s smallest, also dco is smallest. Already in the 
simple diatomic this leads to such a large intra-molecular 
density contribution that the overall p A ™ c is largest. 

These new features also carry over to the correspond¬ 
ing friction coefficients, and in turn to the calculated 
lifetimes. As apparent from Fig. |T] the LDFA-AIM ap¬ 
proach cures the systematic overestimation of lifetimes 
and yields values that are now en par with the FP cal¬ 
culations and experiments. Similar findings have been 
obtained when applying the LDFA-AIM to two further 
systems |47| for which reference lifetimes have been re¬ 
ported from orbital-dependent theories, CN on Pt(lll) 
f5U] and H 2 on Ru(0001) [T7j- With r AIM = 0.9 ps vs. 
r FP = 2.4 ps and r AIM = 210 fs vs. r 1 ' = 215 fs, re¬ 
spectively, LDFA-AIM yields in both cases lifetimes that 
are fully consistent with the reference numbers, yet at 
a fraction of the numerical cost of the orbital-dependent 
theories. 

Adjusting p e mb to also take into account influences 
from atoms other than only the clean metal surface hence 
seems to provide a simple, but effective correction for a 
molecular treatment within the LDFA. The idea underly¬ 
ing Eq. ([3]) is thereby similar to the subtraction of a free 
atom density as suggested in the context of vibrational 
damping of adatoms [44j . However, employing Hirshfeld 
sharing functions [45] offers the advantage that the em¬ 
bedding densities are guaranteed to be within physically 
well defined boundaries 0 < < pscF(R-i)- This 

way, our proposed scheme also preserves the molecular 
dissociation limit by construction: At large bond dis¬ 
tances, the respective friction coefficients smoothly go 
over into friction coefficients virtually identical to the 
ones obtained for independent atoms. This is nicely il¬ 
lustrated by essentially identical vibrational lifetimes ob¬ 
tained within the IAA and AIM for the H 2 on Ru(0001) 
system, where the individual atoms are separated by 
about 2.7 A on the surface. Furthermore, a Hirshfeld 
analysis typically requires only a minute computational 
effort compared to achieving self-consistency for elec¬ 
tronic energies and forces (or compared to the multiple 
self-consistent calculations required to obtain derivatives 
from finite differences in the FP or KT approaches). The 
LDFA-AIM scheme proposed here can thus be easily car¬ 
ried out at every time step of ab initio MD simulations, 
allowing surface motion to be explicitly taken into ac¬ 


count. In this respect, Eq. ([3| defines an embedding 
density not only for adsorbate but also for bulk atoms. 
Friction coefficients derived therefrom could thus, in prin¬ 
ciple, as well be invoked to evaluate energy losses due to 
electron-phonon coupling in the bulk. 

In conclusion, we have shown that the vibrational 
damping of high-frequency adsorbate modes on metal 
surfaces can be added to the list of non-adiabatic phe¬ 
nomena that are reasonably well described by means 
of electronic friction. The complete neglect of intra¬ 
molecular effects in the prevalent LDFA-IAA approx¬ 
imation is thereby likely to underestimate the elec¬ 
tronic energy dissipation. The here presented atoms-in- 
molecules alternative instead accounts for them approxi¬ 
mately through a charge partitioning scheme. As it thus 
effectively treats the molecular electrons as part of the 
metallic substrate, we expect the AIM friction concept to 
generally rather overestimate non-adiabatic energy losses 
and to perform best for chemisorbed adsorbates at close 
distances to the surface. Being a direct descendant of the 
LDFA, our scheme is, of course, also unlikely to overcome 
fundamental limitations that come with its heritage. As 
such it is unlikely to properly capture the strong enhance¬ 
ment of friction coefficients directly at the transition state 
leading to molecular dissociation. However, as important 
as it may be, this actual dissociation event only consti¬ 
tutes a fraction of the dynamics that is relevant to chem¬ 
ical surface reactions. Other important aspects like a vi¬ 
brational pre-excitation or an ensuing hot adatom motion 
may dominate the overall (non-adiabatic) energy dissipa¬ 
tion IS SHI, yet take place over much longer time-scales. 
They can thus only be assessed with a numerically highly 
efficient method like the LDFA. In this respect, our re¬ 
sults further consolidate the trust in LDFA-based results 
for these important long-term events. 
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